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Results of the recent studies of the multiplicity fluctuations and Bose-Einstein 
correlations (BEC) in deep-inelastic scattering (DIS) at large Q'^ are reviewed. 
The measurements were done with the ZEUS detetor at HERA. 

1 Introduction 

Various aspects of multihadron production can be revealed through multiplic- 
ity fluctuations and particle correlations studies. The fluctuations are sensitive 
to soft, particle dynamics and their measurements allow to_test perturbative 
DCD El and local parton hadron duality (LPHD) hypothesis □. In BEC studies 
the correlation function gives information about the shape, size and lifetime 
of the boson emitter source. 

In the fluctuations studies, the method of normalized factorial moments, Fq, 
was used □. The moments of order q were calculated by counting n, the 
number of charged particles in a restricted region of phase space, O: 

Fg{n) = < n < n{n - 1) . . .{n - q + 1) >, (1) 

where < . . . > denotes averaging over the sample and Q is defined either as 
a polar-angle ring around the jet axis or as an upper limit on the particle 
transverse or absolute momentum calculated w.r.t. the jet axis. 

2 Angular multiplicity fluctuations 

The factorial moments were measured as a function, of the ring width 9 in 
the corresponding cone with half opening angle Qq B. Substituting 6 with a 
scaling variable, z = ln(Oo/6')/ ln(iJ0o/A), the experimental results can be 
compared with analytic calculations where E is energy of an outgoing quark 
radiating the gluons. The QCD + LPHD prediction for Fg is: 



In = z (1 - Dq){q 1) ln(i?eo/A), (2) 



where Dq are the Renyi dimensions Q whiclL can be calculated either in a 
fixed or in a running-coupling regime (see u) of the Double Leading Log 
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Figure 1. (Left): Comparison of the factorial moments with QCD+LPHD calculations. The 
curve DLLA (4) includes a correction from the MLLA in the calculation of Dq; (Right): 
Comparison of the factorial moments with different MC models at hadron and parton level. 



Approximation (DLLA). For independent particle production: Dq = 1 and 
Fq{z) = Fq{Q)- Figure l(left) compares the factorial moments for the DIS 
data with the QCD predictions. A significant disagreement with the data 
was found. Figure 1 (right) shows comparison with Monte Carlo predictions 
at hadron and parton levels. All MC models reproduce the trends seen in the 
data. For consistency with the LPHD picture, the parton cascade was cut-off 
at Qq = 0.27 GeV, which is close to A = 0.22 GeV. For higher order moments, 
the parton level is closer to the data than the analytic calculations. 

3 Multiplicity fluctuations in limited momentum space 

According to QCD+LPHD predictions 0, the normalized factorial moments 
are expected to behave as 

when particles are restricted in either the transverse momentum pt < pf"* 
or in spherical momentum p < p'^"*. If p™* Qq, then all Fq ^ 1 and 



ismdSlsd: submitted to World Scientific on February 1. 2008 2 



ZEUS ZEUS 




p" GeV p" GeV p" GsV p" SeV 



Figure 2. (Left): Comparison of the factorial moments calculated as a function of p^"* with 
different MC models at hadron level and parton-level ARIADNE, which represents analytic 
calculations and is consistent with LPHD; (Right): The same for p*^"* factorial moments. 

the multiplicity distribution approaches Poissonian distribution due to the 
coherence effect of soft gluons in a parton cascade. This is not the case 
for soft ghions with spherical momentum cut. Figure 2 shows p™* and p*^"* 
moments El together with Monte Carlo predictions at the hadron and parton 
levels. The parton level is consistent with analytic calculations. For pt < 1 
GeV, Fig. 2 (left), all moments rise in contradiction to the analytic QCD 
predictions. The MC results for hadrons show similar trends to the data. The 
similar behaviour is found for p'^"* factorial moments for small momentum cut 
off, Fig. 2(right). The observed disagreement between QCD+LPHD and the 
measurements is on qualitative level. 

4 Bose-Einstein correlations at large 

Bose-Einstein correlations between pairs of identical bosons can be ex- 
pressed by the two-particle correlation function of the Lorentz-invariant four- 
momentum transfer Q12 — \/^{pi — ^2)^ {pi, P2 are the four- momenta of 
the particles) : R{Qi2) = p{Q 12) / Po{Q 12), where po(Qi2) represents the two- 
particle density without the Bose-Einstein effect. After the corrections for 
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Figure 3. The radius r and strength A of BEC as a function of Q^. 



dynamical correlations and detector effects (Monte Carlo without BEC), the 
correlation function R was fitted to the following expression: 

R{Qi2) - ^(1 + eOi2)(l + Ae-'-''3-), (4) 

where r estimates the size of the two-boson emitter which is taken to be 
of Gaussian shape, A measures the BEC strength, factor 1 + eQi2 takes into 
acount possible long-range momentum correlations and k is the normalization 
constant. 

The BEC were studied to test of the energy dependence of r and A. Figure 
3 shows r and A as a function of four- momentum transfer -squared Q^. No 
dependence was found for r and A. The HI DIS results eI at lower are 
consistent with ZEUS data. 
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